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i <""| We consider top quark pair production in association with a hard jet through next-to- 

i 1 leading order in perturbative QCD. Top quark decays are treated in the narrow width ap- 

' proximation and spin correlations are retained throughout the computation. We include 

hard jet radiation by top quark decay products and explore their importance for basic kine- 
matic distributions at the Tevatron and the LHC. Our results suggest that QCD corrections 
and jet radiation in decays can lead to significant changes in shapes of basic distributions 
and, therefore, need to be included for the description of ttj production. We compare the 
shape of the transverse momentum distribution of a top quark pair recently measured by 
the DO collaboration with the result of our computation and find reasonable agreement. 



I. INTRODUCTION 

• T—i - 

X 

Experiments at the LHC are in the process of accumulating a large data set of top quark pairs 
that will allow detailed studies of various processes that Tevatron experiments either observed with 
relatively low statistics or did not observe at all. Such processes include associated production of 



a tt pair with a jet 



ll], a photon 



3]> t 

tt pair production at very high invariant masses, detailed investigations of associated production 
processes will mark the beginning of the post-Tevatron era in top quark physics. A significant body 
of theoretical work is devoted to improving predictions for tt associated production processes, see 



Refs. 



It is well-known that, once produced, top quarks decay very rapidly. For this reason top quarks 
are observed and studied indirectly through kinematic features of their decay products. Unfor- 
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tunately, this complicates top quark studies by introducing additional uncertainties in kinematic 
reconstructions due to finite resolution on energies and angles of decay products, missing energy 
as well as backgrounds, including combinatorial ones. On the positive side, the rapid decay of top 
quarks enables the description of their decay products in perturbative QCD without the need to 
resort to fragmentation functions and other non-perturbative objects. 

A precise description of hard hadron collisions requires the application of perturbative QCD 
through next-to-leading order (NLO) in the expansion of the strong coupling constant. The com- 
plete NLO QCD description of any process that involves tt production should include QCD cor- 
rections to top quark pair production and to top quark decays. For processes where top quarks are 
produced in association with a photon or a jet, a standard process to study is tiX production with 
X = followed by the top quark decay t — > bW. However, since both photons and jets can be 
radiated in top quark decays, one should also consider tt production followed by "radiative" decays, 
such as t — > bWj and t — > bW^f. The importance of radiation in the decays strongly depends on 
the selection criteria that are used to isolate a particular process and, hence, can not be quantified 
a priori. For example, in a recent measurement of ti'y production by the CDF collaboration Q], 
about half of all signal events come from the process pp — > tt followed by the radiative decay of the 
top quark t — > Wb^f j|. To compare their measurement with theoretical predictions, CDF uses 
a NLO QCD A'-factor for the process pp — > tij computed with stable top quarks. However, since 
about half of their events come from ti production followed by radiative decays of top quarks, it is 
unclear if such a comparison is meaningful. 

In principle, one can get around the problem of separating production and decay stage by 
simply giving up on the approximation that top quarks are produced on-shell and focusing instead 
on the fully realistic final state such as bbW + W~ X with X = "f,j,jj, H, Z. A calculation of pp — > 
bbW + W~ X through a given order in the perturbative expansion in QCD leads to a prediction for 
a final state that includes both "resonant" and "non-resonant" contributions, providing a complete 
description of the process. Without a doubt, this is the best approach possible, provided that it is 
feasible. The feasibility depends on the approximation in perturbative QCD at which the process 
of interest is considered. At leading order, this approach can be pursued for essentially arbitrarily 
complicated process thanks to automated programs such as Madgraph [13] . However, this approach 
becomes very complex already at NLO QCD. For the simplest process pp — > W + W~bb that, among 
many other ways, can occur through the production of a nearly on-shell ti pair, this was recently 
accomplished in Refs. Applications of this approach to more complicated processes are 

difficult to imagine. On the contrary, a sequential treatment of various production and decay 



3 



stages based on the double resonant approximation for t and i can be generalized to processes of 
significant complexity, at least as a matter of principle. This double resonance approximation is 
parametrically controlled by the ratio of the top quark width to its mass Tt/mt ~ 10 -2 and should 
be sufficiently accurate for most observables. In fact, there has been significant progress in using 
this approximation to describe top quark pair production recently. For example, ti pair production 
at NLO QCD in the double resonance approximation, including corrections to top quark decays 
and spin correlations, was computed in Refs. 



16H23I], The number of similar computations for more 



complicated processes is rather 
is associated production of ttj 



imited. The only process for which a full description is available 



where NLO QCD corrections to the production and decays, 



including the radiative one (t — > Wbj), are computed. 

The production of ttj at NLO QCD was first studied in Ref. 
and later in Ref. 
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251 ] for stable top quarks 



261 ] where 



process is described in Refs. 



decays were included at leading order. A different approach to this 
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where ttj production at NLO QCD is combined with a parton 



shower, following the POWHEG procedure 



291 ] . Top quark decays are treated in the parton shower 



app roximation where tt spin correlations are omitted either at leading [27j or at next-to-leading 



281 ] order, and whose correspondence with NLO QCD computations is not clear. 

Fortunately, these approximations are not necessary, since it is possible to treat the complete 
process ttj — > bbW + W~ j in the narrow width approximation where top quark decays, including 
t — > Wbj, are described consistently at NLO QCD and spin correlations are retained throughout 
the entire decay chain. Such a calculation gives a state-of-the-art description of the ttj production 
that, in principle, can be directly compared to experimental results because theoretical predictions 
for a complete and fully realistic final state become available. The goal of the present paper is 



261 ] by including radiation 



therefore to extend the description of pp — > ttj production given in Ref. 
in the decay through next-to-leading order in perturbative QCD. 

The paper is organized as follows. In the next Section, we outline the framework of our calcu- 
lation and discuss technical aspects of the computation which arise because of the need to treat 
radiative corrections to processes with decay kinematics. Phenomenological results for the Tevatron 
and the 7 TeV LHC are presented in Section 3. We conclude in Section 4. 



II. TECHNICAL ASPECTS OF THE CALCULATION 



In this Section, we summarize the technical aspects of the calculation. We begin by describing 
various contributions that we require for the computation. As we pointed out already, the top 
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(a) jet emission in production (b) jet emission in decay 




(c) mixed contribution 



FIG. 1: NLO QCD corrections to top quark pair production and decay in association with a jet. Contribu- 
tions (a) and (b) show jet emission in production and decay, respectively. The symbol "real" indicates that 
one parton is allowed to be unresolved, (c) defines the "mixed" contributions. 

quark is treated in the narrow width approximation. This allows us to organize the computation 
in terms of a production process which includes the hard collision, and the decay process. 
To give a complete list of all necessary contributions for it + jet production calculation, we begin by 
writing the formula for the inclusive cross-section as a convolution of the production cross-section 
a t i and the decay rate Y t 

ddinci = r t ; t 2 ot {da t t +0j + da a+lj + da tl+2j + ...) ® (dr tt - +0j + dr tt - +lj + dr tt - +2j - + ...) . (1) 

Subscripts denote the number of exclusive jets defined according to some jet algorithm. We further 
use the abbreviation dr t j +nj - = X^=o <^t+ij dTj + r n _i)j to summarize the decay rates of top and 
anti-top quark in association with a fixed number of jets. 

We can now expand Eq.(TjQ) assuming that the number of jets that we eventually require is 
equal or larger than one and that the cross-sections and widths for each jet multiplicity scale as 
a tt,nj ~ ^( a s +n ) an d ^t,nj ~ O(a^). Since we are interested in NLO QCD corrections to one-jet 
production, we can disregard all terms that depend on powers of a s higher than four. We obtain 

(2) 

and we re-write this formula in a way that separates various processes that contribute to the 
cross-section 
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(a) 

d^SS = r t ; t 2 ot (dogydi^ + d^drgy + '(d^ + d^ydr^f (3) 
+ d^° (drg y + drgSw) + d ^ dr ^m, + d^ rt dr^ + dagydi^ 

^ V V 

(6) (c) 

We now review different contributions that appear in Eq.Q. The first and second term describe 
ti+j production at leading order followed by leading order decays of the top quark and ti production 
followed by a radiative decay of the top quark, respectively. The third term represents the NLO 
QCD correction to the production process ti + j, where the symbol "real" indicates that one 
parton is allowed to become unresolved. The first term in the second line of Eq.([3]) describes 
leading order production of a top quark pair followed by NLO QCD corrections to the "radiative 
decay" t — > W + b + j . Finally, the last three terms describe mixed contributions where jet 
emission occurs simultaneously in both production and decay stage. Since one of those jets can 
be unresolved, the last two terms are the corresponding virtual corrections needed to provide an 
infra-red finite result. In the remainder of the paper we will refer to contribution (a) and (6) in 
Eq.([3|) as jet radiation in the production and jet radiation in the decay, respectively. The last part 
(c) we call the mixed contribution. The corresponding topologies are depicted in Fig. [TJ 

Let us now describe how NLO QCD corrections to jet radiation in the production processes 
pp — > ti and pp — > tij are treated. We note that - when production processes are considered 
at next-to- leading order - the decay processes are included at leading order, consistent with the 
expansion in a s . However, these leading order decays are different processes: in the former case, 
we consider the radiative decay t — > Wbg, since an additional jet is required in the final state. In 
the latter case, top quarks decay into the Wb final state since the jet is created in the production 
stage. The NLO QCD results for the production processes are available; they are described in 



Refs. 
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261 ] including an efficient way of implementing the decays of top quarks while retaining 



all spin correlations. We note that one-loop QCD corrections to — > qqtt, — > ggtt, — > qqttg and 
— > gggtt amplitudes that we require are calculated using generalized D-dimensional unitarity 



3C- 



32l | . The real emission corrections are obtained following the Catani-Seymour dipole subtraction 



formalism 



33l | and its extension to massive particles in Ref. [3^ . To improve the efficiency of the 



computation, we follow Ref. [37f] and use a-parameters to restrict subtraction terms to singular 
phase-space regions. The relevant dipoles with a-parameters are found in Refs. 
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36] 



1 We note that, in the case of a semi-leptonic top quark decay, also the W-boson is allowed to radiate an additional 
hard jet at NLO QCD. We include this contribution in our computation as well. 



ti 



The second part required in the calculation involves leading order production processes pp — > ti 
and pp —> tij followed by top quark decays at next-to-leading order. In the former case the 
NLO QCD corrections to radiative decays t — > Wbj and W —> qqg are required; in the latter 
case t —> Wb and W —¥ qq need to be computed through NLO QCD. Radiative corrections to 
t — > Wb and W — > qq are known; our implementation follows the description in Ref. 22|. We 



do not repeat it here and focus, instead, on the NLO QCD corrections to the "radiative decay" 
t — > bWj. Since this is a sufficiently low- multiplicity process, we compute the virtual corrections 
using Passarino-Veltman reduction of tensor integrals 
Ref. 39[]. For the calculation of the real corrections we need to consider various decay processes, 



381 ] . The scalar integrals are taken from 



such as t — > (W qq') bgg, t — > (W — > qq'gg)b and t — > (W — > qq'g)bg etc. The real emission 



subtraction terms are again constructed using the dipole formalism of Catani and Seymour |33l |. 
However, we note that its application to decay processes requires clarification. Catani and Seymour 
constructed subtraction terms - the dipoles - that satisfy two criteria: 1) they remove infra-red and 
collinear singularities when subtracted from scattering amplitudes and 2) they can be integrated 



331 ] . it is shown how to satisfy 



analytically over the unresolved phase-space. In the original paper 
these conditions for two colliding massless partons. Since decay kinematics differ from production 
kinematics, some of the Catani-Seymour dipoles need to be modified if we deal with decays of 
color-charged particles. 

Recall that within the Catani-Seymour dipole formalism, dipoles are constructed by taking 
different partons to be "emittors" and "spectators", in addition to soft or collinear partons that 
are actually "emitted". The dipoles depend on "flavors" (quarks, gluons) of "emitted" and "emit- 
tors" and on whether "emittors" and "spectators" are in the initial or in the final state. The 
corresponding dipoles are referred to as final-final, final-initial, initial-initial and initial-final. 

However, only a limited number of these dipoles is needed for the decay processes in general. 
First, it is obvious that there are no initial-initial dipoles since there is just one particle in the 



331 ] and the phase-space re-mapping 



initial state. Final-final dipoles can be borrowed from Ref. 
therein. Initial-final dipoles can be omitted since real radiation by a massive initial state particle is 
only singular in soft kinematics. This contribution can be absorbed into final-initial dipoles which 
are the only dipoles for decay kinematics that need to be constructed. 

The complete list of dipoles that we need for the process t — > Wbg±g2 are 

^9i32,6) ' D b gi) g 2 , ^bg 2 ,gi ^b g v ^692 and ^9192 ■ Tlie first three dipoles are of the final-final type 
whereas the last three dipoles are the missing final-initial dipoles. We will discuss their construc- 
tion in the following. We need to distinguish two types of final-initial dipoles which correspond to 
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the splitting q — > qg and g ^ gg with a top quark in the initial state being the spectator. 

We begin our discussion with the gluon-quark dipole. It can be extracted from Ref. [jcj]. To 
this end, we consider the process t — > Wbg\g2 and imagine that gluon g\ and the (massless) b- 
quark become unresolved. The top quark in the initial state is the spectator. We combine the 
momenta of the VF-boson and the gluon 52 into a new momentum p\y = pw +p g2 and introduce a 
variable r 2 = p\^jm\. The remaining momenta - whose scalar products lead to soft and collinear 
singularities - are parametrized using two variables z and y 



PbP 9l 



(1 



r) 2 y, 



PtPgi 



ml 



(l-r 2 )(l-z). 



With this parametrization, the final-initial gluon-quark dipole reads 



4-7ra s /x 



2c 



1 



PbP 9l 



1 



1 - z - ye(l - z) 



mi 



(PtPgi, 



(4) 



(5) 



where e = (4 — d)/2 is the parameter of dimensional regularization, d is the number of space-time 
dimensions and A, A' are quark helicity labels. We note that Eq.([5]) gives the dipole in conventional 



dimensional regularization (CDR) scheme; if four-dimensional helicity (FDH) scheme 
the term proportional to e in Eq.(j5]) should be dropped. 



411 ] is used, 



In Ref. [9| we have integrated the dipole in Eq . ([5|) over the restricted unresolved phase-space [37j] , 



drawing extensively from the results of Ref. 



401 ]. We reproduce this result here for completeness. 



We consider the integration of the dipole in Eq.(j5]) over the unresolved restricted phase-space 



[dg] [i 



l- a- z)0{y - ay max )] D* l6 

-y)~ e [1 - 6(1 - a - z)9{y - ay 



l/max 

x J dyy~ e (y 




1 

M J dz (r 2 + z(l 




(6) 



ixj] ^gib- 



where 



(1 +r) 2 z(l - z) 
z + r z (l — z) 

We find the following result in CDR 



TV 



j [dg] D t gib [l-e(l-a-z)e(y 
a s (4vr/i 2 ) e 



16vr 2 



"2/max)] 



-m 



2-2e 
t 



(47T) 



r(i-e) Vi 



1 + r 



2c 



(7) 



2vr mf r(l - e) 



8\\>, 



1 1/5 , , ,\ 27 1/1 
- + - - - 2 In(n) + - + - 3-- + 7)lnH 



''1 



1 ^-TT-^ 

— + 2Li 2 (n) - — - 51n(n) +21n 2 (n) 
lr\ o 



-2 In 2 a 



1 — n 



7 , « 2 \ 2(1 -a)r 2 , 

4a H In a + — — In , 

2 2i n VI -an 



(8) 
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with r\ = 1 — r 2 . 

It remains to construct the gluon-gluon dipole of the final- initial type for decay kinematics. 
In variance with the gluon-quark dipole just considered, the gluon-gluon dipole contains non- 
trivial spin correlations. We will use the parametrization of the unresolved phase-space that we 
just discussed with an obvious modification of the momentum pw', for the gluon-gluon dipole, it is 
given by pw = Pw+Pb- To derive the gluon-gluon dipole, we consider the limit of the — > tbg\g2W 
amplitude squared when two gluons become collinear. The result reads 



\M\ 2 ^MIPHM\ (9) 



where 



P99 ^ 



-^( T f ? + - r i ]-2(l-e)«l-0-^ 



(10) 



is the spin-dependent splitting function. In Eq. lfTOj) , £ and are defined as 

(l-£)(2pn)' Pg2 qp 1 £(2pn)' 
where the light-like vector p defines the collinear direction and another light-like vector is 
auxiliary. We can now use the relations between gluon momenta 

k 1 ~ <V = 6& - (1 - 0P&, 2p 9lP52 = -^Ly, (12) 

to write 



1 fiu 



( £ ,l _ , ft \ a v a » 

- s '-{— ( + —) +{l - f - ) M 

To construct the dipoles, we split this expression into two terms 

jd99 

1 ILL/ 



(13) 



^ ~ D 1 ' 2 + D 2 ' 1 (14) 



where 



ni,2 = j WOL I x _ c " " 

^ 2 Ws2 \ a-e) + 2 ( Ma )i 1 j 

and -D^'i is given by Eq. (ll5p with £ — > 1 — £. We would like to rewrite this equation in such a 
way that the integration over the unresolved phase-space becomes straightforward. To this end, 
we express Eq. (|15p in terms of the variables z and y and momentum of the top quark pt and pw- 
Because 

(PtP 9l ) _ 1-g Mfi v 
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we can identify £ with the variable z in Eq.Q. It remains to modify the spin-correlation part of 
Eq. (ll5p and write it in appropriate variables. We note that such modifications can be arbitrary 
provided that the original form of the spin-correlation part of the dipole is recovered in the limit 
when p gi and p 92 become collinear. We do that by writing 

a ^^ = - ptPi2+PtP^ \ au _ (17) 

V PtPl2 J 

In Eq. (|17|) . the momentum p±2 is the light-like vector given by p\2 = Pt — A-pw, where A is 
the Lorentz transformation constructed explicitly in Ref. {^(J. The reduced matrix element that 
describes the decay process t — > W+b+g is then evaluated for p t , Apw, and pu, where Ap\y is then 
split into the W momentum and the 6-quark momentum. We note that the projection operator 
introduced in Eq. (|17p ensures that ir^ is transverse to p±2- As we show below, this feature simplifies 
the integration over the unresolved phase-space considerably. It is straightforward to check that 
in the collinear (y —> 0) limit, tt^ — > an. Hence, to construct a suitable dipole, we can simply 
substitute ir^ for in Eq. f)15|) . Note also that we are allowed to multiply the spin-correlation part 
in Eq. (|15p by an arbitrary function f(y,z) provided that it is free of singularities and that it is 
normalized in such a way that /(0, z) = 1. We choose this function to be 

tl \ 4 {PtPw) 2 ~ r 2 mf 

f(y, z) = — — - — ^2 — ' ( 18 ) 

to simplify the calculation of the integrated dipole with a-dependence. As the very last step, we 
add one more term to the dipole, to account for soft singularities that appear when a gluon is 
emitted from the top quark in the initial state. We are finally in the position to write down the 
gg final-initial dipole. In the CDR scheme, the result reads 

1 



D )j,u,t _ a 2e 

^PgiPgv 

The various quantities that appear in Eq. (fT9j) are 
1 



au , ^_ _ m t 2p giPg2 \ (1 - e )vr^ 

\i- z 4 (p tPgi y) + 2 PgiPg2 ny ' z \ 



(19) 



(20) 



((ptPu)^ - J^(pi 2 a)) , a» = [(PtPffaKi - (PtP gi )Pg 2 

PtPn m tv- ~ r ) 

771? . 9 . , , 777,?, N n 777? (1 — r 2 ) 

PtPgi = " r 2 )(l - z), p giPg2 = - r) 2 y, p t p 12 = '-, 

with r 2 = (pw +pb) 2 /m 2 . 

To integrate the dipole in Eq. (119j) over the unresolved phase-space, we make use of the results 
presented in Ref. [^(J. It is straightforward to integrate the part proportional to the metric tensor. 
Integration of the spin-correlation part is more involved but it can be simplified because vector 7r^ 
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is orthogonal to the light-like vector p^ 2 - This allows us to write 



(2P1P2) 



Ai 



y. - 



PtPl2 



+ A2P12P12 



(21) 



where (•••) 3/ , 2 denotes the integration over y and z as in Eq.(jSJ). The term proportional to A2 
can be dropped since it gets Lorentz-contracted with the product of on-shell matrix elements that 
vanish when contracted with p\2- Hence, we only need to compute A±, which we easily obtain by 
contracting the left hand side of the above formula with the metric tensor. By the same argument, 
once A\ is obtained, we can drop terms proportional to p^ 2 in tensorial structure that is multiplied 
by Ai in Eq. (|2ip . Therefore, we can write the result of the integration of Dg"'g 2 over unresolved 
phase-space as proportional to the metric tensor. 

We now present the result for the integrated final-initial gg dipole in the CDR scheme for decay 
kinematics, including its full a-dependence. The integrated dipole reads 

(47T/X 2 ) 6 



[dg] D££ [1-9(1- a- z)9(y - ay m3X )] 



2vr mf r(l - e] 



9^u 



1 17- 12 logri 
2e* + 12e 



^ , 2 

log a 

12 B 



(1 - a) (23 - a + 2a 2 ) 



17 1 

- -j- log n 

6 



r log r 



12 

2, 



log q + log ri 



6rf 

+ 12an(r 2 + rf) + r 2 (2 + n(-l + lln)) 



6a 3 (l - ri)(-2 + n) - 3^(1 - ri)(-6 + 5n) 



+ 



'1 — a)r 2 log(l — ar±] 
4rf 



(-2a 2 (l - n)(-2 + n) + a(-2 + (5 - 3n)n) - 2 + n + r 2 - 4rf ) + Li 3 (n) 



240r^(l - an) 



8a 9 rf 



6a 8 rf (2 - 7n) - a 7 r 3 (20 - 68ri + 115r 2 ) 

+ a 6 r 2 (-40 + 130n - 165r 2 + 216r?) - a 5 ri(120 - 360n + 410r 2 - 234r 3 + 305r£ 
+ a 4 (240 - 180ri - 510r 2 + 650r 3 - 195rf + 278rf) 
- a 3 (600 - 1140ri + 280r 2 + 460r 2 - 92rf + 97r?) 
+ a 2 (360 - lUOn + 900r 2 + 50r 3 - 63rf - 40rf) 

+ 10a (12 + 6ri - 36r 2 + 10r 3 + 8r? + 91rf) + 10r 2 (4r 2 - 91r 2 ) 



(22) 



The integrated dipole given in Eq. (|22p is the final ingredient we need to treat the real emission 
contributions to radiative decays of top quarks. 
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III. PHENOMENOLOGICAL RESULTS 



In this Section we present phenomenological results for the Tevatron (y/s = 1.96 TeV) and the 



LHC (y/s = 7 TeV). We choose m t = 172 GeV for the top quark mass and m\ 



/ 



80.419 GeV 



for the VF-boson mass. We employ MSTW2008 parton distribution functions [42J and use the 
corresponding values of a s at leading and next-to-leading order. The couplings of the W-boson 
to fermions are obtained from the Fermi constant G F = 1.16639 • 10 _5 GeV~ 2 . Since we work in 
the narrow width approximation, our results are inversely proportional to the top quark and the 
W-boson widths, a ~ l~y 2 r^?. These decay widths are evaluated at leading and next-to- leading 
order in the strong coupling constant, for LO and NLO cross-sections, respectively. For reference, 
we give the results for the widths 

= 1.4653 GeV, Tf LO = 1.3375 GeV, 

(23) 

= 2.0481 GeV, T^ LO = 2.1195 GeV. 

The shown NLO results for the widths are computed with the renormalization scale fx = mt- We 
note that the use of NLO expressions for the widths increases the NLO cross-sections by about ten 
percent. 

We begin with the discussion of the Tevatron results. We consider tt production in the lepton 
+ jets channel so that our leading order cross-section contains five jets. The lepton transverse 
momentum and the missing energy in the event are required to satisfy p±i > 20 GeV and E™ 1SS > 
20 GeV. Jets are defined according to the fcj_-jet algorithm [43J with AR = 0.5. The jet transverse 
momenta are required to be larger than pj_ j > 20 GeV. Both leptons and jets must be central 
\yi\ < 2, \yj\ < 2. To better discriminate against the background, we require an additional cut on 
the transverse energy in the event H± = ^2jP_i_,j + p±, e + E™ 1SS > 220 GeV. We present results 
below for a single lepton generation. Hadronic decays of VF-bosons to first two quark generations 
are included and the CKM matrix is set to the identity matrix. 

The cross-sections for pp — > bW + {e + v e ) bW~(jj) + j production at the Tevatron at leading and 
next-to-leading order in perturbative QCD, subject to the above cuts, read 

ctlo = 75.291^-2 fb) f7NLQ = 78 . 9= 5.6 ^ (24) 

In Eq . (|24|) . the central value refers to renormalization and factorization scales set to fj, = m t and 
the upper (lower) value to fi = m t /2 and \i = 2mt, respectively. We observe a dramatic reduction 
in dependence on unphysical scales if NLO QCD corrections are included. 
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20 40 60 80 100 120 140 160 180 200 220 240 260 
PrQet) [GeV] 



FIG. 2: Fractions of events when the leading (non-6) jet at the Tevatron comes from tij production, the 
decay t — !> Wbj or mixed processes, as a function of jet transverse momentum. Note the sign of the mixed 
contribution and the cancellation between decay and mixed mechanisms at high transverse momentum. 
Rcnormalization and factorization scales are set to fi = m t . 

It is interesting to understand how jet radiation in the production and jet radiation in the decay 
contribute to cross-sections shown in Eq. (j24p . To answer this question, we present separate cross- 
sections for production and decay processes as well as mixed contributions, as defined in Eq.Q. 
For factorization and renormalization scales set to fi = mt we find 

a LO = 46.33 (Pr) + 28.96 (Dec) = 75.29 fb, 

(25) 

ct nlo = 47.7 (Pr) + 36.7 (Dec) - 5.5 (Mix) = 78.9 fb. 
This result is interesting because it shows that, with our choice of selection criteria, in only sixty 
percent of all events that contain a ti pair and a jet, the jet can be associated with the production 
process; in the remaining forty percent of events, jets come from top quark decays. These fractions 
are stable against NLO QCD corrections, but the reason for that stability is peculiar. Indeed, it 
follows from Eq. ([25p that the NLO QCD corrections to the production process are relatively small 
[K = 1.03) while QCD corrections to the decay process are quite large {K = 1.37). There is, 
however, a significant negative contribution from the "mixed" corrections. As described around 
Eq.([3|), this contribution arises from single jet emission in the production convoluted with single 
jet emission in the decay and the corresponding virtual corrections. Because of this cancellation 
between decay and mixed contributions, a relatively small correction to jet radiation in top quark 
decays remains. Thus, an estimate of the NLO cross-section that employs the exact leading- 
order cross-section as in Eq. (|24p and the fT-factor for the production process K = 1.03 gives 
o"lo X K = 77.54, which is in good agreement with the full NLO result (/i = mt) in Eqs. (124)25p . 
However, this cancellation seems accidental to us. In spite of the proximity of the two numbers 
for the tij production at the Tevatron, we were unable to come up with a convincing and general 
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FIG. 3: Distributions of the lepton transverse momentum, the lepton rapidity, the transverse momentum 
and the rapidity of the hardest jet for ttj production at the Tevatron at leading and next-to-leading order in 
perturbative QCD. The bands correspond to the variation of rcnormalization and factorization scales in the 
interval rrit/2 < fx < 2m t . Results with hard jet emission in the production stage only followed by leading 
order decays t — > W + b are compared to full NLO results in lower panes. 

argument that ensures that .fT-factors for the production and decay processes are always similar. In 
fact, the importance of mixed and decay contributions strongly depends on the kinematic variables. 
For illustration we show production, decay and mixed contributions as the function of the transverse 
momentum of the leading non-6 jet in Fig. [2j At low p 1 ? < 60 GeV, jet radiation in top quark 
decays is the largest (~ 60%) contribution to the cross section. As expected, at larger , the 
jet is predominantly emitted in the tt production. The mixed contribution is positive at small jet 
momenta but changes sign at moderate p 1 ? and cancels the contribution due to jet radiation in 
decay at large . The situation appears to be quite complex and observable-dependent. We can 
therefore anticipate - and we will see this explicitly in the context of the LHC discussion - that 
calculations without accounting for jet radiation in the decays of top quarks can lead to misleading 
results. 

Various kinematic distributions at the Tevatron are shown in Figs. [3] and 01 For all kinematic 
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FIG. 4: Distributions of the transverse momentum and the rapidity of the 5th hardest jet, the transverse 
energy Ht and the transverse momentum of the ti pair for ttj production at the Tevatron at leading and 
next-to-leading order in perturbative QCD. The bands correspond to the variation of renormalization and 
factorization scales in the interval mt/2 < /i < 2mt- Results with hard jet emission in the production stage 
only followed by leading order decays t — > W + b are compared to full NLO results in lower panes. 



distributions we find a significantly reduced dependence on the choice of the factorization and the 
renormalization scales as well as shape changes in kinematic tails of some distributions. The impact 
of QCD radiation in top quark decays is illustrated in the lower panes of each plot, where ratios 
of full NLO cross-section and the NLO ttj production cross-section followed by the leading order 
decays of top quarks are shown. In general, these plots confirm the expectation that QCD radiation 
in top quark decays mostly affects spectra at low transverse momenta. But there are interesting 
exceptions where the impact of radiation in the decay is more pronounced. In particular, we 
find fairly uniform enhancement of transverse momenta and rapidity distributions of the charged 
lepton as well as the rapidity of the hardest jet (Fig. [3]). The decay contribution to the rapidity 
distribution of a lepton is asymmetric; it appears to be more important at large positive rapidities. 
However, the full NLO distribution does not show significant asymmetry in lepton rapidity. 

In Fig. H] we show distributions of the transverse momentum and rapidity of the 5th hardest 
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FIG. 5: Comparison_ of DO measurement with theoretical NLO QCD prediction. The data points are 
obtained from Ref. 



44[ after background subtraction. The bands correspond to the variation of renormal- 
ization and factorization scales in the interval mj/2 < fx < 2m t . The experimental distribution and the 
(j, = n%t theoretical distribution are normalized in such a way that their integrals equal to one. 

jet, the total transverse energy in the event H± and the transverse momentum of the ti pair. All 
these distributions receive non-uniform enhancements from jet radiation in top quark decays. In 
particular, H±_ and pj_(5th jet) distributions are strongly enhanced at low values of H± and p±, 
where relatively soft radiation in top quark decays dominates. Also, the rapidity distribution of 
the 5th hardest jet receives strong enhancement at central rapidities which is a consequence of the 
fact that top quark decay products are produced mostly at small rapidities. We note that similar 



shape changes were recent 
approximation in Ref. 



y observed in the context of studying pp — > ttj within the parton shower 



n 



271 ] . Note, however, that the cross-section computed in Ref. [27[ seems 
closer to the contribution that we identify as "jet radiation in production". While - as we just saw 
- such a result underestimates the cross-section, it is probably consistent with the fact that decays 
in Ref. are treated in the parton shower approximation which by construction conserves the 
overall probability and does not change normalization. 

We also consider the distribution in the transverse momentum of the ti pair in Fig. [H This 
kinematic distribution is particularly interesting because recent results by the DO collaboration [44 ] 



show a disagreement between predictions of MC@NLO [45J] and data at low transverse momenta. 
Since we deal with top quark decay products rather than with stable top quarks, we need to define 
what is meant by the ti transverse momentum. To this end, we imagine that the reconstruction 
proceeds by finding two non-6 jets whose invariant mass is closest to Myy and then combining 
the transverse momenta of these two jets, two 6-jets, the lepton transverse momentum and the 
missing transverse momentum, to obtain the transverse momentum of the tt pair. We find that 
the transverse momentum distribution of the ti pair is affected by the radiation in the decay 
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FIG. 6: Fractions of events when the leading (non-6) jet at the 7 TeV LHC comes from ttj production, the 
decay t — > Wbj or mixed processes, as a function of jet transverse momentum. Note the sign of the mixed 
contribution and the cancellation between decay and mixed mechanisms at high transverse momentum. 
Renormalization and factorization scales are set to fi = m t . 

non- uniformly - the decay contributions are more important for small values of p± (tt) . 



44| , we combine the pp — > ttj 
22j. In the pp -> tt 



To further compare the results of our computation with DO data 
calculation described above with a pp — >■ tt computation at NLO QCD 
computation we impose a jet veto prohibiting additional jets with the transverse momentum larger 
than 20 GeV, for consistency with the current tij computation. We present our results 2 and the 



DO data 



441 ] in Fig. [5j The normalization of the fi = mt NLO computation is chosen such that the 



integrals of the two distributions agree. In spite of the significant theoretical uncertainty in the 
lowest bin, it appears that our calculation can well describe the shape of the p±(tt) distribution 
observed by the DO collaboration. The lower pane in Fig. [5] shows that inclusion of NLO QCD 
corrections to pp — > ttj and pp — > tt is crucial for achieving the agreement. 



We continue with the discussion of tij production at the y/s = 7 TeV LHC. We imagine that W- 
bosons from both t and t decays decay leptonically. For definiteness, we assume that the top quark 
decays to a positron and the antitop quark decays to an electron. All generic input parameters that 
we employ in the calculation were already described at the beginning of Section III. Specific to the 
LHC case, we require at least three jets, defined by the anti-fej_ jet algorithm Q with AR = 0.4. 
All jets have a minimum transverse momentum p±j > 25 GeV and central rapidities \yj\ < 2.5. 
Similarly, leptons need to satisfy p± j > 25 GeV and \yi\ < 2.5, and the missing energy in the event 



2 We note that the kinematic cuts on the final state particles that we use are similar but not identical to the ones 
used by DO collaboration. 
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FIG. 7: Distributions of the lepton transverse momentum, the lepton rapidity, the transverse momentum 
and the rapidity of the hardest jet for tij production at the LHC (7 TeV) at leading and next-to-leading 
order in perturbativc QCD. The bands correspond to the variation of renormalization and factorization 
scales in the interval mt/2 < fj, < 2mt- Results with hard jet emission in the production stage only followed 
by leading order decays t — > W + b are compared to full NLO results in lower panes. 

^miss y Q e y yy e fi nc [ following results for leading and next-to- leading order cross-sections 

<7 LO = 350.3+^-0 ^ aNLO = 2g8 -46 ^ (2g) 

In Eq. (l26p . the central value refers to renormalization and factorization scales set to [i = mt and 
the upper (lower) value to fi = m t /2 and fi = 2mt, respectively. 

In case of the LHC, the interplay between radiation in the production and radiation in the decay 
is very different from the Tevatron. Since top quark pairs at the LHC are mostly produced in gluon 
annihilation and the collision energy is high, radiation in the production strongly dominates over 
radiation in the decay. We find (/i = mj) 

ct lo = 316.9 (Pr) + 33.4 (Dec) = 350.3 fb, 

(27) 

CTNLo = 323 (Pr) + 40.5 (Dec) - 75.5 (Mix) = 288 fb. 
The three NLO contributions are shown in Fig. [6j as a function of the transverse momentum of 
the leading non-6 jet. The radiation in the decay becomes less and less important as the process 
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FIG. 8: Distributions of the transverse energy Ht, the transverse momentum of the top quark pair, the 
di-lepton invariant mass and the relative azimuthal angle between the leptons for ttj production at the LHC 
(7 TeV) at leading and next-to-lcading order in perturbativc QCD. The bands correspond to the variation 
of renormalization and factorization scales in the interval mt/2 < fx < 2mt- Results with hard jet emission 
in the production stage only followed by leading order decays t — > W + b are compared to full NLO results 
in lower panes. 



becomes harder, but the negative mixed contribution appears to be significant also at high pj_. 
Although radiation in the decay at the LHC is less important than at the Tevatron, it is peculiar 
that "mixed" contributions are large and negative. 

We point out that this may cause misleading results, if the full (production and decay) lead- 
ing order cross-section and the next-to-leading .fT-factor for the production process only are used 
to estimate the full NLO cross-section. The iiT-factor (/i = mt) for the production process is 
323 fb/316.9 fb ~ 1.02, so the naive estimate of the NLO cross-section is 1.02 x cjlo ~ 357 fb, 
which is about twenty percent higher than the correct NLO value given in Eq. (l27p . We emphasize 
that the "mixed" contribution to ttj production is a NLO QCD effect, so unless NLO effects are 
properly incorporated into computations of associated production of unstable particles, it is unclear 
to what extent various predictions for cross-sections can be trusted. 

In Fig. [7] and [8] we show various kinematic distributions for the LHC. The importance of QCD 
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radiation in decays for various observables can be seen from the lower panes. We find that for 
the LHC, the impact of the QCD radiation in the decay is modest; the variable that seems to 
be most affected is H± at small values of the transverse energy. For kinematics distributions in 
dilepton invariant mass or in the relative azimuthal angle of the two leptons, there is a uniform 
reduction, almost independent of mi+i- and <Pi+i-- Finally, we note that given the discrepancy 
between MC@NLO prediction for the transverse momentum of the tt pair and the DO data 44 ) . it is 



important to measure this distribution at the LHC. Thanks to a much higher energy and luminosity, 
the LHC should be able to probe a much broader distribution in p±(tt), including regions where 
fixed order QCD computations are directly applicable. We show the p±(tt) distribution in Fig. [8] 
and find that this distribution receives important modifications due to radiation in the decay. 



IV. CONCLUSIONS 



In this paper, we discussed the computation of NLO QCD corrections to the production of a 
ti pair in association with a hard jet at hadron colliders. While NLO QCD corrections to this 
process have been considered in the literature several times already, in this article for the first 
time, QCD radiative corrections to top quarks decays are studied, including the possibility that 
the jet is emitted in the decay stage. The results reported in this paper lead to a complete and 
fully consistent treatment of top quark pair production and decay in association with a jet at 
next-to-leading order in perturbative QCD. 

While at leading order there is a clear separation into production and decay stages, at next-to- 
leading order there appears a new contribution where one parton is emitted in the production and 
the other parton in the decay. Since this "mixed" contribution must be supplemented by virtual 
corrections to ensure infra-red safety, we find that it can be negative. This leads to interesting 
effects that, to the best of our knowledge, have not been discussed in the literature before. In 
particular, it is far from clear that a widely used procedure of estimating NLO QCD cross-sections 
by computing leading order cross-sections with decays and re-scaling them by if-factors obtained 
from calculations that ignore radiation of jets in the decay is valid. In fact, we find that this 
procedure accidentally gives an accurate estimate of the NLO cross-section for ttj production 
at the Tevatron but similarly overestimates the NLO QCD cross-section at the LHC by twenty 
percent. The absence of clear pattern suggests that it is best to include QCD radiative corrections 
to decays of unstable particles into theoretical predictions for hard scattering processes. 

Jet radiation in the decays can have significant impact on kinematic distributions. One such 
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case is the H± distribution at the Tevatron which exhibits significant distortion due to radiation 
in the final state. While the situation at the LHC is less dramatic, even there certain distributions 
are systematically distorted at the ten to twenty percent level. 

We also compare the shape of the transverse momentum distribution of a top quark pair recently 
measured by the DO collaboration with the result of our computation. We combine exclusive 
pp ti and inclusive pp — > tij computations at NLO QCD to describe the transverse momentum 
distribution of ti pair and find reasonable agreement with the results obtained by DO collaboration 
in Ref. [4| . 



Recent progress in NLO computations was driven by the idea that perturbative QCD can 
describe hard scattering well, pushing theorists towards providing realistic descriptions of com- 
plicated hard processes which can be directly compared to experimental data. Clearly, in the 
case of heavy short-lived particles such as top quarks, this implies that NLO QCD computations 
should be applied to their decay, including all spin correlations. All of this can be done in a 
rather straightforward way in the narrow width approximation which provides a parametric 
framework for such studies. We have demonstrated how this framework can be used to describe 
the production of ti pairs in association with a jet at hadron colliders. We look forward to further 
applying this framework for the description of both Standard Model and New Physics processes 
at the LHC. 
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